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Abstract. We report the spectroscopic confirmation of a compact Lyo- emitting galaxy at z = 5.721. A F0RS2 spectrum of the 
source shows a strong asymmetric line with a flux of 5x10"" erg s"' cm"^, making it one of the brightest Lya emitting galaxies 
at this redshift, and a line-of-sight velocity dispersion of 400 km s"'. We also have a tentative detection of a second, narrower 
component that is redshifted by 400km s"' with respect to the main peak. A F0RS2 image shows that the source is compact, 
with a FWHM of Cf .'S, which corresponds to 3.2 kpc at this redshift**. This source is a brighter example of J1236.8-I-6215 
jPawson et all20o3) . another Lyo- emitting galaxy at z ~ 5.2. 

Key words, galaxies: high-redshift - galaxies: evolution - galaxies: starburst - galaxies: individual (Jl 14334.98-014433.9) 



1. Introduction 

Wide-field imaging surveys with specially selected narrow- 
band filters are an effective means of discovering high redshift 
{z > 5) Lya emitting galaxies (see |Hu et al. 2004^ and refer- 
ences therein). Spectra of these galaxies are dominated by a sin- 
gle, asymmetric emission line. One of the strongest arguments 
for associating this line with Lyman-a (Ly a) is the asymmetry 
in the line profile (e.g. IStern et aObOOOl) . which can only be 
detected if the spectral resolution is high enough (R > 2000). 

The profile of the Lya line is the end result of emission 
from Hn regions and resonant scattering by Hi. The bulk of 
the Lya emission comes from the recombination of hydrogen 
that has been ionised by UV flux of massive stars. Part of the 
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* Based on observations taken at the Cerro La Silla (ESO programs 
67.A-0063, 68.A-0363, 69.A-0314 and MPG time) and Cerro Paranal 
Observatories (ESO program 272.A-5029). 

** Throughout this Letter, a cosmology with Hq = 70 km s"' Mpc"' , 
£1m = 0.3 and CIa = 0.7 is assumed. Magnitudes are on the AB-system. 



ionisation may be due to shocks iIb land-Hawthorn & Nulser] 
I2OO4I) or from an A GN, although in a study of Lya emitting 
galaxies at z ~ 4.5. IWang et al.l (|2004) found no evidence of 
AGN activity. The shape of the line profile is sensitive to the 
geometry, density and kinematics of both the Hii gas, where the 
Une is produced, and the Hi gas, where it is scattered jAhn et all 
l2003llSantos et alj2004l) . Dust can also play a role. Hence, the 
Lya emission line can be viewed as a tool, albeit a rather blunt 
one, that might be used to constrain the spatial and kinematic 
distribution of the hydrogen gas in these distant galaxies. 

In this letter, we present the (J.'S resolution seeing-limited 
imaging and R ~ 3600 spectroscopy of a Lya emitting galaxy 
that was selected from Lya emitting candidates in the WFILAS 
catalog (Westra et al. 2005, in prep.). 

2. WFILAS and Candidate Selection 

WFILAS is a survey for bright Lya emitting galaxies at z ~ 
5.7. The selection strategy is similar to the successf ul strategies 
employed in other Lya surveys at these redshifts (lAiiki et al.l 
f2003t iRhoads et al...2003. .Hu et al...2004) . WFILAS covers a 
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Fig. 1. From left to right, thumbnails of the confirmed Lyff emitter SI 1_13368 at redshift 5.721. Each thumbnail covers a 24" x 
24" region with a pixel scale of 0'.'238 pix"' . The emitter is only detected in the M^is and A^sn images. 



larger volume (~1.2xlO*Mpc^'), a larger area (~1 n°) and has 
a brighter detection limit (2cr limit magnitude ~24.0-24.5) than 
these surveys. Hence the candidate Lya emitting galaxies in the 
WFILAS catalogue will, on average, be brighter than the Lya 
emitting galaxies in these other catalogs. 

The survey used the Wide Field Imager (WFI) on the 
ESO/MPI 2.2m telescope at the Cerro La Silla Observatory 
and targeted three fields. The WFI consists of a mosaic of 
eight (4 X 2) 2k x 4k CCDs arranged to give a field of view 
of 34' X 33' with a pixel scale of 0'.'238 pixel"'. Images were 
taken with the standard broad-band B and R filters, and four 
narrower filters - an intermediate-band {AA = 22 A) filter at 
815 nm (Mgis) and three custom-made narrow-band (AA = 7 A) 
filters with central wavelengths at 810nm (A^bio), 817 nm (A^sn) 
and 824 nm (A^824)- They lie in a spectral region where the emis- 
sivity of the night-sky is relatively low, which improves the sen- 
sitivity to Lyo' emission. 

Candidate Lya galaxies are those that appear in one of the 
narrow band filters, but are undetected in the broad band fil- 
ters. Given the relatively low signal-to-noise ratios of the can- 
didates, we also require a detection in the intermediate band 
filter This limits the number of spurious candidates. We refer 
the interested reader to Westra et al (2005, in prep.) for a de- 
scription of the observations, reduction and candidate selection. 
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Fig. 2. A 30" X 30" region around the confirmed Lya galaxy 
S11.13368 atz = 5.721 from the pre-image taken with FORS2. 
The image has a pixel scale of 0'.'252pix"'. The object is unre- 
solved in this image. The seeing at the time of these observa- 
tions was ~ 0'.'5. The exposure time was 3600 seconds. 



3. Confirmed Lya Emitter at z = 5.721 

In a pilot study to test the effectiveness of the selection strategy 
one of the brighter candidates (Jl 14334.98-014433.9, hereafter 
SI 1_13368) was observed with FORS2 on Yepun (UT4) at the 
Cerro Paranal Observatory (Figure^. 

A pre-image with an intermediate-band filter (13 nm) cen- 
tered at 815 nm was taken with FORS2 on 2004 February 16th 
in which SI 1_13368 clearly was detected. Figure|2lshows a 30" 
region around SI 1.13368. The FWHM of stars in this field are 
0'.'5 and Sll_13368 is unresolved. 

Three 1200 sec exposures were taken on 2004 March 18th 
with FORS2 using the 1028z grism and a 1" slit. Frames 
were bias-subtracted and flatfielded and were then combined 
with suitable pixel rejection to remove cosmic rays. The 2D- 
spectrum (without subtracting the sky lines) and the extracted 
sky-subtracted spectrum are shown in Figure|3fl, b), where one 
can clearly see a single emission line with a broad red wing. No 
continuum is detected, implying a 2cr upper-limit for the con- 
tinuum of 7x10"^° erg s"' cm"^ A"' over the rest frame wave- 
length range 1220 to 1230 A. 



IStern et alJ tOOd) have reviewed the different ways in 
which high redshift Lya can be verified and suggest that 
line asymmetry is the surest way. Our spectral resolution 
(R ~ 3600) is high enough to securely confirm this asymme- 
try. This resolving power is also more than adequate to rule out 
[On] /1/13726,3728 at z = 1.19, since we do not resolve the 
line into the close doublet. The separation of this doublet at 
z = 1.19 is 6.1 A and therefore easily resolvable. Similarly, we 
can rule out the possibility that the line is [Oiii] /1/14959,5007 
at z = 0.63, due to the absence of the accompanying line in that 
doublet. We can also rule out Ha at z = 0.25. If the line was Ha, 
then we should have either detected [N n] and/or some flux in 
the /?-band due to the contributions of H/3, [Oiii] .1.14959,5007 
and the continuum (e.g.^niazev et al. 2004)- In Figure^ one 
can see that there is no detection in the iJ-band image. Given 
the absence of all of these potential neighbouring features, and 
also the clear asymmetry of the line, we identify it as Lya emis- 
sion at z = 5.721. 

The integrated line flux of the line derived from the spec- 
trum is 5x10"'^ ergs"' cm"-, making it one of the bright- 
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Fig. 3. A fit of both one- and two-component models to the 
profile of the Lya line, (a) The two dimensional spectrum cen- 
tered on the Lya line. In this unrectified spectrum, the night 
sky emission hnes are not removed. The pixel scale for this 
image is 0'.'252pix"' in the spatial direction and 0.86 Apix"' 
in the dispersion direction, (b) Observed Lya line (histograms) 
with the three best-fitting models. For clarity, the models are 
off'set by 3 erg s"' cm"^ A"' . The two-component models con- 
sist of a two components: a broad, truncated Gaussian and a 
narrower redshifted Gaussian. The one-component model con- 
sists of only a broad, truncated Gaussian. See Tabled for the 
parameters of the 'broad' (solid), 'spiky' (dotted) and 'single 
component' (dashed) models, (c) Same model line profiles as 
in (b) but before convolution with the instrument profile, (d) 
Observed data minus model fit (as plotted in (b)) residuals, 
demonstrating a random scatter about the zero flux line. Also 
shown (histograms) is the Icr-error spectrum from the observed 
data, which includes both sky- and Poisson noise. Note that the 
red peak is not N v. At this redshift it would appear around 
8334 A. 

est Lya emitting sources at this redsh ift (cf. lAiiki et al ] |2003[ 
iRhoads et al.1l2003llMaier et al.ll2003l) . The 2cr lower limit on 
the rest frame equivalent width is ~100A. At z - 5.721, this 
translates to a Lya luminosity of 1.8x10'*'' ergs"' suggesting 
an apparent star- formation rate of 16 Moyr using the con- 
version rate from lAiiki et alJ (l2003h. 

Fo llowing earlier works ( IPawson et alJ l2002i iHu et all 
I2OO4I) . we fit both two- and single-component models to 
the Lyo' line. The two-component fit consists of a truncated 
Gaussian with complete absorption bluewards of Lya line cen- 



ter an d a redshifted Gaussian that is not truncated (e.g. lHu et alJ 
I2OO4I) . The one-component fit consists solely of a truncated 
Gaussian (Figurel^J^.c). Since the seeing was narrower than the 
width of the slit, we convolve the model with a Gaussian that 
has a FWHM of 2.3 A. At 8175 A this corresponds to a resolu- 
tion of R ~ 3600. We use the Levenberg-Marquardt nonlinear 
least-squares algorithm to find the best fit. 

Two different two-component models fit the data with a 
similar reduced values of ~ 1.3. We refer to these two mod- 
els as the 'broad' model and the 'spiky' model. In the 'broad' 
model, the redshifted component is broader and weaker and the 
central component is narrower and stronger in comparison to 
the 'spiky' model. After convolving with the instrumental pro- 
file both fits have similar residuals. Both models have a central 
peak at a wavelength corresponding to Lya at z = 5.721. The 
redshifted component is clearly detected in both models and 
lies ~ +400 km s"' away from the central peak. Given the sim- 
ilar reduced we cannot favour one model over the other. 

The single-component model has a broader main peak, 
which is slightly bluer This model does not fit the profile as 
well as the two-component cases, particularly in the region of 
the red peak (xl - 2.2). 

Table[2summarises the different model components. 



Component 
(1) 


Ac /peak FWHM 

(2) (3) (4) (5) 


Av 

(6) 


(7) 


Main peak 
Red peak 


'Broad' model 

8173.1 7.0 11.1 

8184.2 1.9 2.2 


408 
81 


+ 11.1 


406 


Main peak 
Red peak 


'Spiky' model 
8173.1 6.9 11.4 
8184.1 4.3 0.6 


419 
24 


+ 11.0 


405 


Main peak 


'Single component' model 
8172.2 6.8 14.0 514 







Table 1. Model fit parameters as described in Section|3]and in- 
dicated in Figure|3l Notes: (1) component of the fit, (2) central 
wavelength of the fitted component in A, (3) peak flux density 
in 10""* erg s ' cm ^ A"', (4) and (5) FWHM of fufl Gaussian 
of the profile in A and kms"', respectively, (6) and (7) line-of- 
sight outflow velocity in A and km s"', respectively. 



The integrated line fluxes for the two-component models 
are very similar. If one were to include the flux that was missing 
from the blue side of the truncated Gaussian, the total a line flux 
is 8.3x10"'^ erg s"' cm"-. This corresponds to a sta r-formation 
rate o f ~ 27 M© yr"', using the conversion rate from lAiiki et alJ 
(•2001. The width on the blueward side of the profile is solely 
due to instrumental broadening. 

4. Discussion 

We have presented a medium resolution spectrum of a bright 
Lya emitting galaxy at z = 5.721. The spectrum consists of a 
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single emission line and no continuum. The line shows a dis- 
tinct asymmetry, which undoubtedly confirms it as Lya. We 
model the line with two components: a one-sided Gaussian and 
a narrower, redshifted component. The profile of the blue side 
of the line is entirely defined by the instrument profile. 

Generally, the second component is less frequently ob- 
served, although it is possible that it has been missed in the 
spectra of other Lya emitting galaxies. Most of these spec- 
tra were taken at lower resolution and are considerably nois- 
ier. A second peak in the Lya line is a clear signature of 
an expanding sh ell of neutral hydrogen jPawson et alJl2002l 
lAhn et alj2()03l) . The strength and shape of the secondary peak 
depends on the kinematics and the quantity of neutral hydrogen 
in the expanding shell and the am ount and distribution of dust 
throughout the galaxy (lAhnl2004 . 

Sll_13368 appears to be a brighter and more distant ex - 
ample of J1235.8+6215 at z = 5.190 dOawson et alJl2002l) . 
The line profile is strongly asymmetric in both objects, and 
both suggest a second redshifted component. Both objects are 
also very compact. However, there are some noteworthy differ- 
ences. The Lya line in SI 1_13368 is considerably broader, and 
the redshifted component is a lot narrower, even in our 'broad' 
model. 

The intrinsic Lya profile is heavily modified by the sur- 
rounding gas and the frac tion of the line that is fina lly observed 
is very model dependent ilSantosl2nn4lAhnl2nn4 . In general, 
it is only a fraction of the intrinsic flux. Hence, star formation 
rates that are estimated from the observed Lya flux directly, as 
they are done in this paper, could drastically underestimate the 
true star formation rate. Similarly, the centroid of the observed 
profile is also model dependent. This directly leads to an uncer- 
tainty in the redshift of about 0.01, if no other lines are visible, 
which is usually the case for such high redshift galaxies. 

Sll_13368, like J1235.8+6215, is very compact. With a 
projected size of ~ 3 kpc or less, it is comparable to the size 
of the star forming regions in local starbursting galaxies; how- 
ever, the star formation rate is much higher. Not all Lye emit- 
ting galaxies at z ~ 5.7 are as compact. In Figure 5, we plot 
apparent size of the star forming region versus the inferred 
star formation rate for a sample of local starbursts and distant 
galaxies. The emission line region in LAE J1044-0130 occurs 
over a region that is an order of magnitude larger than emis- 
sion line regions in SI 1_13368 and J1235. 8+6215 even though 
the inferred star formation rate is significantly less. Given that 
the projected star formation rate per unit area in Sll_13368 
far exceeds 0.1 Moyr"' kpc"^, it is likely that a hot, enriched 
starburst-driven gas is outflowing into the halo of Sll_13368, 
facilitating the enrichment of the halo and the e scape of Lyman 
continuum photons (iHeckman et alji2000i iTenorio-Tagle et aD 
Il999h. 
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Fig. 4. Size of the star-forming reg ion versus the star-form ation 
rate. The SFR of two galaxies from lHeckman et al](ll990l) is in- 
dicated by a bar, as there are two measurements of the SFR. The 
arrows represent lower-limits to the SFR and upper-limits to 
the size of the major axis, (due to seeing limited observations). 
Error-bars are not included. The pl us signs represent thre e well- 
known local starbursting galaxies ( lHeckmanetalJl990l) . 
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